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Abstract 

In the present paper we continue our reconsideration about the foundations for a thermostatistical description of the called 
Hamiltonian nonextensive systems. After reviewing the selfsimilarity concept and the necessary conditions for the ensemble 
equivalence, we introduce the reparametrization invariance of the microcanonical description as an internal symmetry associated 
with the dynamical origin of this ensemble. Possibility of developing a geometrical formulation of the thermodynamic formalism 
based on this symmetry is discussed, with a consequent revision about the classification of phase-transitions based on the 
concavity of the Boltzmann entropy. The relevance of such conceptions are analyzed by considering the called Antonov 
isothermal model. 
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I. INTRODUCTION 

As alsewhere discussed, the study of the non extensive 
systems constitutes an interesting and fascinating chal- 
lenger for the developments of the Thermodynamics and 
Statistical Mechanics 0,lllllllllillllll|lSlil|ll 
Ha. Il4| . An important feature of this kind of systems 
is the existence of long-range correlations, whose origin 
delays on the existence of long-range interactions among 
the constituents particles or the small (or mesoscopic) 
nature of the system. 

We address in our previous paper [13 a reconsidera- 
tion of the foundations for a thermo-statistical descrip- 
tion of the special case of the nonextensive Hamiltonian 
systems. The main conclusions derived from this work 
are the following: (1) The universality of the microscopic 
mechanism of chaoticity supports the generic applicabil- 
ity of a thermostatistical description with microcanonical 
basis for the nonintegrable many-body Hamiltonian sys- 
tems [ia, lis HM HM I2O1 im J even for those nonextensive 
systems with long-range interactions 22, 23, 24, 23, (2) 
Extensive conditions of the traditional systems are not 
longer applicable for the nonextensive systems, but such 
properties could be generalized by considering the scal- 
ing selfsimilarity concept 15] which is considered in order 
to establish an adequate thermo-statistical description in 
this context. 

We will continue in the present paper our analysis of 
the foundations for a thermo-statistical description of 
the nonextensive Hamiltonian systems. Our aim now is 
to show how the scaling selfsimilarity concept could be 
used in order to perform an appropriate thermodynamic 
formalism in the systems with long-range interactions. 



As already discussed in our previous work, most of the 
Hamiltonian systems with a practical interest exhibit ex- 
ponential selfsimilarity scaling laws [I^- This is the rea- 
son why the present discussion will be focussed in this 
special class of the nonextensive Hamiltonian systems. 
Interested reader could also review a similar approach to 
the one developed in this work in ref. 26] for the case 
of the potential selfsimilarity scaling laws, which seem 
to be associated with the popular Tsallis nonextensive 
statistics Wi- 



ll. REVIEWING SELFSIMILARITY 

Selfsimilarity is just a general symmetry of the macro- 
scopic description exhibited under the scaling transfor- 
mations of the system size in a times, N —^ N (a) — 
aN, which manifests itself in every macroscopic behav- 
ior of the system, even in the macroscopic dynamics ll5| . 
From the mathematical viewpoint, selfsimilarity appears 
as a multiplicative uniparametric group of transforma- 
tions Tq, TaTfs = Tap, which acts on those fundamental 
macroscopic observables / determining the microcanon- 
ical macroscopic state of a given system, leading in this 
way to a scaling transformation of the number of micro- 
scopic configurations W = W (!) (microcanonical acces- 
sible phase-space volume) whose form does not depend 
on the macroscopic state and can be represented as fol- 
lows: 



/ (a) =TaI=>W [a] = y {ay-^ (W)) 



(1) 
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where the function y (x) determines what we refer as the 
selfsimilarity scaling law of the system. According to the 
framework developed in the previous work 15], the self- 
similarity scaling law y (x) determines the generalized 
Boltzmann Principle (microcanonical entropy) and the 



specific statistical entropic compatible with the selfsimi- 
latity of a given system: 






(2) 



where L {x) = y^^ {x). Such definitions should be con- 
sidered as the starting point in the generalization of the 
traditional thermodynamic formalism. 

Two remarkable cases of selfsimilarity scaling laws are 
the exponential selfsimilarity: 



y{x) = exp(a;) ^ Sbg = -^Pfclnpfc, 



(3) 



which leads to the usual Boltzmann-Gibbs statistics, and 
the potential selfsimilarity: 



y (x) = expg {x) ^ Sq ^ ~^pI lUqPk, 



(4) 



associated with the Tsallis nonextensive statistics |27|. 
where In, x = exp~^ (x) = (a;^~' — l) / (1 — q), is the q- 
generalized logarithmic function. Among them, the most 
extended in practical applications is the exponential self- 
similarity scaling laws © . Extensivity of traditional sys- 
tems is just a special case of this kind of scaling selfsim- 
ilarity, whose scaling selfsimilarity transformations obey 
the well-kown form: 



N (a) ^ aN 
E (a) = aE 
V (a) = aV 



W{a) =exp(alnTy). 



(5) 



While a generic Hamiltonian system of the form: 



with the total potential energy V {q), i,j = l,2,...n 
and q = (^q^, q^ , . . . g"), exhibits an exponential growing 
of the microcanonical accessible phase-space volume W 
with the increasing of the system degrees of freedom n 
[l5|. supporting in this way an exponential growing of 
W during the scaling operation n —f n (a) = an, the 
selfsimilarity scaling transformations T^ acting on the 
fundamental macroscopic observables (such as the total 
energy E or the systems volume V) are far to look like the 
extensive case (|S)| in the long-range interacting systems. 
A very simple example of system with nontrivial scal- 
ing transformations is found for the selfgravitating gas 
of identical nonrelativistic particles with a finite incom- 
pressible volume V which interact among them by means 
of the Newtonian gravity. Obviously, the particles den- 
sity p in this system obeys the inequallity p < po = 1/w. 
For very low energies and large densities where p ~ po 
the scaling transformations are approximately given by: 



while for very high energies and low densities where the 
particles size can be dismissed since p ^ po, such scaling 
transformations can be expressed asyntotically as follows: 

N (a) = aN, E (a) ~ a^E, V (a) ~ a'W, (8) 



results which have been obtained in refs.j23, |22| and it 
will be shown again in the section Ivl 

Generally speaking, the determination of the specific 
selfsimilarity scaling transformations T^ could be a very 
difficult task because of there is no methodology to es- 
tablish them in a given application outside the extensive 
context. As already suggested in the introductory sec- 
tion of our previous paper [13 , systems with long-range 
interactions look like the extensive systems in the neigh- 
borhood of the second-order phase transitions due to the 
long-range correlations existing in this region. This is 
the reason why we think that the well-known Renormal- 
ization Group methods 30] should play a crucial role in 
the determination of the scaling selfsimilarity of a given 
nonextensive system. 



III. THERMODYNAMIC FORMALISM 

As elsewhere discussed, the thermodynamic formalism 
is based on the equivalence in the thermodynamic limit 
N —^ 00 between the microcanonical ensemble ujm and 
the canonical ensemble lubg |31| : 



1 



ujm 



nii,N) 



5 I-Ia 



1 



N I , ^BG 



Z{f3,N) 



exp 



{-PIn) 



(9) 

derived from the consideration of the Maximal Entropy 
Principle |32l | by using the extensive entropy Q, where 
In are all those macroscopic observables (integrals of mo- 
tion) determining the macroscopic state. Obviously, the 
thermodynamic limit N ^ 00 is equivalent to carry out 
a scaling transformation Ta on the observables / with 
a ^- 00. 

The ensemble equivalence can be shown by considering 
the Laplace transformation which relates their partition 
functions: 



Z(/3,7V)= / e-f^'n{I,N)dI. 



(10) 



Introducing the Boltzmann entropy Sb {I, N) = In W, 
where W = flSI^, and the Planck thermodynamic po- 
tential V (/?) = —\nZ (/3) , the Laplace transformation 
could be rewritten as follows: 



e^v{-V{P,N)) 



exp I 



[PI- 



Sb{I,N)])^. 



(11) 



N (a) = aN, E (a) ~ a^E, V (a) ~ aV, 



^ Here 51 is a small volume constant associated with the coarsed 
(7) grained of the phase-space. 



We immediately recognize in the argument of the expo- 
nential function the well-known Legendre transformation 
between the thermodynamic potentials: 



P™ (/3™, iV) = /3™/ - Sb (/, N) , with /3,„ = 



OSb 
dl 



(12) 



where the subindex m means that such quantities are 
associated with the microcanonical description. 

It can be shown by using the steepest descend method 
that the equivalence between Vm {(3m, N) and V {P,N) 
takes place in the thermodynamic limit when the follow- 
ing conditions are applicable: 

i) The macroscopic observables / behave in a 
extensive-like way during the selfsimilarity scaling 
transformation T^: 



Tal ^al ^ Sb {T^I, aN) = uSb {I, N) . 



(13) 



ii) There is only one stationary point i^ = I/N satis- 
fying the relation: 



(3 = 



dsB (is) 
dl '■ 



(14) 



where sb (i) — j^Sb {iN, N) is the entropy per par- 
ticle. 

The condition i) is only trivially satisfied by the ex- 
tensive systems when the microcanonical description is 
performed by using the called extensive observables, i.e. 
total energy, electric charge, etc. However, extensivity 
does not take place in systems with long-range interac- 
tions. In this context extensivity should be substituted 
by the scaling selfsimilarity, but the extensive-like char- 
acter H13|l of the macroscopic observables / does not hold 
generally speaking. 

The condition ii) demands the concavity of the Boltz- 
mann entropy, property which can not be everywhere 
satisfied. The existence of several stationary points leads 
to a catastrophe in the Legendre transformation which 
is related with the occurrence of first-order phase tran- 
sitions at a given value of (3. Such anomaly appears in 
those regions where at least one eigenvalue of the Hessian 
of the Boltzmann entropy: 



Kri 



d^SB 
QinQim ' 



(15) 



is positive. As already shown by Gross, every anomaly in 
the Hessian of the Boltzmann entropy can be related with 
the existence of phase transitions, extending in this way 
the identification of such critical phenomena is systems 
outside the thermodynamic limit 



reparametrization invariance of the microcanonical de- 
scription. ^What is such reparametrization invariance? 
As already commented, the microcanonical ensemble is 
just a dynamical ensemble, and therefore, every symme- 
try of the microscopic dynamics is also relevant at the 
macroscopic level in the microcanonical description. The 
reparametrization invariance is just a symmetry with a 
dynamical origin. 



A. Mathematical definition and Physical meaning 

Let us consider now all those integrals of motion 
/ {X) = {/I {X) ,P{X),... /" {X)} determining the 
microcanonical thermostatistical description of a given 
nonintegrable many-body Hamiltonian system. Each 
admissible value of the above set integrals of motion 
/ = {/^, J^, . . . /"} determine certain subset Sp of the 
phase-space X : 

X eSp={X eX\l''{X) = I'', k:^ 1,2,... n}. (16) 

in which is enclose the system microscopic dynamics in a 
given microcanonical macroscopic state. 

The totality of these admissible "points" {/} defines 
a subset TZj of the n-dimensional Euclidean space 7?.", 
while the totality of the corresponding phase-space sub- 
sets {Sp} defines certain partition 3 {X) of the phase- 
space X: 



'^(X)^ ISpCX 



\^Sp — X ; SpC] Sq — 



(17) 



The above definitions H16(l and (|17|l allow the existence of 
certain bijectivc map ipj between the elements of 3 (X) 
(subsets Sp C X) and the elements oiTZi (points / G 7?.") 

^i : 5 (A") ^ 7^/ = {ySp e 5 (X) 31 eUi C TC 
\iX eSp^ I'' {X) = /^ fc = 1, 2, . . . n} . (18) 

The bijective character of this map provides to the parti- 
tion 3 (X) the same topological features of certain subset 
TZi of the n-dimensional Euclidian space TZ". Now on we 
refer the partition 3 (X) as the abstract space of the inte- 
grals of motions, and it will be denoted simply by 5. We 
say that the map ipi defines a n-dimensional Euclidean 
coordinate representation TZj of the abstract space 3. 

The coordinate representation TZj of 3 is not unique. 
Let TZi and TZ^ be two diffeomorphic subsets of the n- 
dimensional Euclidean space TZ", that is, there exist a 
diffeomorphic map: 



IV. REPARAMETRIZATION INVARIANCE 

We will show in this section that the conditions i) 
and ii) could be satisfied by considering the called 



ip : TZi ^ TZ^ 



yi eTZi 3(fe TZ^ 



i-ll^l^o 



(19) 
among their elements. The terms "diffeomorphic map" 
means that the n-dimensional Euclidean subsets TZi and 



^Z^p posses the same topological and diffeomorphic (differ- 
ential) structure. The map '^|}^p = ilJioip~^ defines a new 
bijective application among the elements of the space 9 
and the points ip € TZ^p, ip, 



■ip, yip 



-*■ TZip^ and consequently, 
Tpip is also another n-dimensional Euclidean coordinate 
representation of 3. The coordinate representations TZi 
and TZp are equivalent because of they represent the same 
abstract space 9. 

Generally speaking, a reparametrization change is car- 
ried out when we move from the Euclidean representa- 
tion TZi towards an equivalent Euclidean representation 
Ti-ip. The totality of the diffeomorphic maps (|19|l consti- 
tutes certain group of transformations called the group 
of diffeomorphisms of the space 9, DiffC^), which is 
the maximal symmetry that a geometrical theory could 
exhibit. 

Every reparametrization change H19() also induces a 
reparametrization of the integrals of motions ip => px '■ 



ipx : i (X) ^ piX) ^ [p' (^i {X)) , . . .p" (^ 



iX))}, 
(20) 

Since / (X) are integrals of motions, every p'^ {I (X)) e 
p {X) will be also a integral of motion. Due to the bijec- 
tive character of the reparametrization change p : TZi — > 
TZip, the new set of integrals of motion p{X) generates 
the same partition 3 (A") of the phase-space. This is the 
reason why the sets / {X) and p {X) can be considered 
as equivalent representations of the relevant integrals of 
motion determining the microcanonical description. 

Since the totality of the phase-space points belonging 
to a given subset 5p G 3 (X) represents the same micro- 
canonical macroscopic state, and the diffeomorphic and 
topological structure of the abstract space 9 does not 
depend on the reparametrization changes, it is not dif- 
ficult to understand that the microcanonical description 
is reparametrization invariant. 

Theorem 1 The microcanonical ensemble ujm 0) is in- 
variant under every reparametrization change p> : 7?./ 
-^Tlp. 

Proof. Such invariance is straightforward derived from 
the following identity of the Dirac delta function: 



8U{I)-ip[h 



N = 



dp 



dl 



S(I-I. 



N 



(21) 



which leads to the following transformation rule for the 
partition function Q: 



nip,;N) 



dp 
'dl 



niI;N), 



(22) 



and consequently: 
1 



S{p- p>n) = 



1 



n{p>,N) '"^ ^ ' n{i,N) 

UJM (v5, N) = UM (/, N) 



S(I-I. 



N 



(23) 



We have considered here the implicit dependence of in- 
tegrals of motion on the system size N. ■ 

A corollary of the Theorem ^ is that the Physics 
derived from the microcanonical description is 
reparametrization invariant since the expectation 
values of every macroscopic observable O derived from 
the microcanonical distribution function ojm exhibits 
this kind of symmetry: 

0= f OcjAidX ^0{p,N) = (/, N) , (24) 

relation which is straightforwardly followed by using the 
identity if^ . 

The reparametrization invariance does not introduce 
anything new in the macroscopic description except the 
possibility of describing the macroscopic state by using 
any representation TZp equivalent to the original repre- 
sentation TZi , a situation analogue to describe the phys- 
ical space TZ^ by using a Cartesian coordinates (x, y, z) 
or spherical coordinates {r,0,p>). 

Notice that the diffeomorphic reparametrization 
changes demand the analyticity of the integrals of mo- 
tions, which is a necessary condition for a particular inte- 
gral of motion /'^ be relevant in the macroscopic descrip- 
tion. On the other hand, the integrals of motions / are 
derived from their commutativity with the Hamiltonian 
H. Since H is one of those integrals of motion determin- 
ing the macroscopic state (total energy), the Hamiltonian 
loses its identity during the reparametrization changes 

(|20|l . The only way to preserve the condition <H,I> — 

during the reparametrization changes is by demanding 
the commutativity among all those integrals of motion 

determining the macroscopic state I I^^r 1 = 0: 



{^^v^^■}-E 



kl 



dp}^ dtp^ 
dl'' dP 



{/^/'} 



(25) 



Interestingly, this is the necessary condition for the si- 
multaneous determination of physical observables in the 
Quantum Mechanics, and therefore, reparametrization 
invariance is also consistent in this level. 



B. Some geometrical aspects 

Reparametrization invariance is a symmetry of the mi- 
crocanonical description which allows us to perform a 
geometrical description in this framework, a covariant 
formulation of the Therniostatistics. Geometrical formu- 
lations of the Therniostatistics are very well-known in 
Physics. Some recently proposed geometrical formula- 
tions are given in refs. Q and [SSLIS^ . 

A very important question is to determine the transfor- 
mation rules of the physical observables and the thermo- 
dynamical function under the reparametrization changes. 



Now on the Einstein summation convention will be as- 
sumed. As already shown, the microcanonical distri- 
bution function and the physical observables behave as 
scalar functions during the reparametrization changes: 

COM (if, N) - LUM (/, N) , Oi^,N)^0 {I, N) . (26) 

The transformation rule microcanonical partition func- 
tion Q corresponds with the transformation rule of a 
scalar tensorial density: 



n{^,N)^ 



dip 



dl 



n{i,N). 



(27) 



The first derivatives of a given physical observable O be- 
have as the components of a covariant vector: 



do 



dl'' do 



gjk 



dip"^ dip"^ dl^ 



dip-- 



-Vk- 



(28) 



The transformation rule of the second-order tensors is 
given by: 



Qjk gjl 

dip"^ dip"" 



A 



kl- 



(29) 



However, such covariant second-order tensor cannot be 
obtained by differentiation without the introduction of 
a Riemannian structure, that is, without introducing a 
Riemannian metric. 

The microcanonical partition function fi allows us the 
introduction of the invariant measure d/i ~ Qdl. Tak- 
ing into account such invariant measure, the integration 
of a given physical observable O defined on the space 
3 (microcanonical ensemble) is directly related with the 
phase-space integration: 



S S(A') 



O (X) dX, 



(30) 



where S is a subset of the space 3, and S {X) is its 
corresponding image in the phase-space. 

Strictly speaking, the Boltzmann entropy Sb = In M^ 
is not an scalar function. The microcanonical accessi- 
ble volume W is just the invariant measure of a small 
neighborhood E of the interest point obtained after by 
performing a coarsed grained partition of the space 3: 



W^fi{E)= J 



= dfi. 



(31) 



In applications W is taken approximately by M^ ~ ildl, 
where the volume SI of the coarsed grain subset S is re- 
garded as a very small constant. The small constant SI 
becomes unimportant when N is very large, that is, with 
the imposition of the thermodynamic limit, and conse- 
quently, the Boltzmann entropy could be considered as a 
scalar function in this limit. 



C. Consequences of this dynamical symmetry on 
the thermodynamic formalism 

The microcanonical ensemble provides a full character- 
ization of the macroscopic properties of a given isolate 
Hamiltonian system in thermodynamic equilibrium. The 
imposition of the thermodynamic limit iV — > oo leads 
under certain conditions to the equivalence of the micro- 
canonical description with certain simplified character- 
ization: the macroscopic description performed by the 
canonical ensemble. As already discussed, such canoni- 
cal ensemble is determined from the relevant statistical 
entropy derived from the specific scaling selfsimilarity 
exhibited by the interest system, i.e.: the exponential 
selfsimilarity is associated with the Shannon-Boltzmann- 
Gibbs extensive entropy lj2l , which naturally leads to the 
called Boltzmann-Gibbs distributions 0. 

It is very easy to verify that the canonical distribution 
function Q does not exhibit the original reparametriza- 
tion invariance of the microcanonical ensemble. This 
means that the Physics in the canonical ensemble de- 
pends on the coordinate representation TZj of the ab- 
stract space 3 of the relevant integrals of motion used 
for perfoming the macroscopic characterization of a given 
system. Therefore, the thermodynamic formalism could 
not satisfy this kind of symmetry, and consequently, the 
ordering information derived from its consideration is not 
reparametrization invariant. 

As already discussed in section lllll the equivalence be- 
tween the microcanonical and the canonical descriptions 
demands the fulfillment the conditions i) and ii), and 
therefore, such ensemble equivalence could not be per- 
formed by using an arbitrary representation. However, 
the existence of the reparametrization invariance of the 
microcanonical descriptions plays a crucial role in arriv- 
ing to a well-defined thermodynamic formalism. 

The extensive-like behavior of the integrals of motion 
under the scaling transformations is simply satisfied by 
chosing an appropriate coordinate representation TZ^ of 
the space 3 exhibiting this kind of behavior under the 
scaling transformations. There are some trivial exam- 
ples in which this demand is very easy to accomplish. 
Considering the example illustrated at the end of section 
im the total energy of the selfgravitating gas of identical 
point particles interacting by means of the Newtonian 
gravity grows with the scaling following a as power law 
|28| . and therefore, an appropriate representation with an 
extensive-like behavior could he p — Np{E/Eo), where 
Eq = GM"^ /R is the characteristic constant energy which 
evidently follows the same as power law, being ip {x) a 
bijective function. 

The reader may notice that such procedure differs from 
the called Kac argument 35j because of the total energy 
is not arbitrarily forced to be extensive, since its scal- 
ing behavior is determined from the system selfsimilarity. 
We only demand the extensivity in the coordinate repre- 
sentation of the relevant integrals of motion in order to 
ensure the ensemble equivalence in the thermodynamic 



limit. The main difficulty in chosing an appropriate rep- 
resentation of the relevant integrals of motion is that the 
specific selfsimilarity scaling transformation Ta acting on 
the space 3 is a priori unknown and have to be deter- 
mined. 

As already shown, the transition from the microcanoni- 
cal towards the canonical description with the imposition 
of the thermodynamic limit and the consequent establish- 
ment of the system selfsimilarity leads to a breakdown of 
the original reparametrization symmetry. However, ^is 
the reparametrization invariance completely broken dur- 
ing this transition? The answer is no. This fact is very 
easy to understand by considering the above example 
about the selfgravitating gas: anyone bijective function 
(fi = Nip (E/Eq) could be used with the purpose of de- 
scribing the macroscopic characterization, and therefore, 
there is a complete set of admissible coordinate represen- 
tations C — {TZip} exhibiting an extensive-like behavior 
during the scaling transformations. 

Let TZi and TZ^p be two coordinate representations be- 
longing to the set of admissible representations C and 
ip = tpipOipy : TZj — > TZip, the reparametrization change 
between them. Since the scaling transformations for 
every coordinate in both representations are given by 
T^{I) — al^ and T'^ [tp) — a(p™, the reparametriza- 
tion change ip (/) is a homogeneous map satisfying the 
relations: 



r 



{al) = a<^™ (/) ^ ^" (/) = /'^^Cr^ 



dl^ 



(32) 



for m — 1, 2, ... 71. 

The complete set of reparametrization changes acting 
in the set of admissible representations C constitutes the 
Group of Homogeneous transformations Ti, which is a 
subgroup of the Group of diffeormorphisms Diff (S) al- 
ready commented. The symmetry associated with the 
transformations of the Homogeneous group Ti. is just the 
residual symmetry of the original reparametrization in- 
variance of the microcanonical description. 

Theorem 2 The thermodynamic potential Vm obtained 
from the Legendre transformation fjji)] is Ti- invariant: 

r^ = l'^^^SsiI)^^-^^-Ss{^). (33) 



dlf" 



dip-" 



Proof. The demostration is straightforwardly followed 
by performing a variable change in the partial derivative 
of the Boltzmann entropy: 



^,dSB ^^,fdip"^dSB 



dl^ 



\ dl^ dip-" 



^,dip"^\ osb 



dl'' J dip-" 



(34) 



and considering the homogeneous character of the 
reparametrization change 1)32(1 : 



I -TTTT- = P 



dl^ 



dip-" 



(35) 



We have shown in this way that the "microcanoni- 
cal" thermodynamic functions Sb and Vrm as well as 
their first derivatives exhibit a covariant behavior un- 
der the reparametrizations changes of the Homogeneous 
group Ti. However, it does not mean that the thermo- 
dynamic formalism exhibits this kind of symmetry be- 
cause of the conditions ii) is necessary to satisfy also in 
order to ensure the ensemble equivalence. The fulfill- 
ment of this condition could be perturbed because of the 
reparametrization changes can affect the concavity of the 
Boltzmann entropy. Let us to illustrate this by consider- 
ing the following example. 

Let ip he a map between two seminfinite Euclidean 
lines, ip : TZ^ — > 7?.+ , which in the representation TZx is 
given by the concave function ip {x) — ^/x (where x > 
0). Let us now consider a reparametrization change "0 : 
TZx -^ TZy given hy y = x^ (which is evidently bijective). 
The map ip in the new representation TZy is given by the 
function ip {y) = y^ (where y > 0), which clearly is a 
convex function. 

The above trivial example teaches us that the concav- 
ity of the Boltzmann entropy depends on the representa- 
tion 7^/ e C used to perform the canonical description, 
and therefore, ensemble equivalence (or inequivalence) 
depends on the coordinate representation of the space 
5. This is a very interesting point since the ensemble 
inequivalence is a signature of the first-order phase tran- 
sitions. 

According to the Microcanonical Thcrmostatistics of 
Gross ^], the Hessian of the Boltzmann entropy H15() 
gives a complete characterization of the phase-transitions 
in the microcanonical ensemble. However, in spite of the 
aparent similarity, the components of the Hessian are not 
the components of a covariant tensor of second-order be- 
cause of they obey the transformation rule: 



d'^Si 



dl^ dP d^S 



d^I^ dS 



dip-"^dip-" dip"" dp-" dfldP dip-""dip"- dl^ ' 



(36) 



while the correct transformation rule is given by 
()29|l . This means that such classification of phase- 
transitions based on the concavity of the entropy is not 
reparametrization invariant. Therefore, the coordinate 
representation of the space 3 should be specified in order 
to precise the existence of first-order phase transitions. 
How such ambiguity in the first-order phase-transition 
existence could be physically acceptable? 

The question is that the coordinate representation TZi 
of the space 3 is determined from the external constrains 
imposed to the system. This idea is very easy to under- 
stand analysing the case of the extensive systems. Ordi- 
narily, canonical ensemble iobg =^ Z^^ {f3) eicp {— P E) is 
experimentally implemented by putting the interest sys- 
tem in a thermal contact with a heat bath. This experi- 
mental arrangement not only fix the system temperature 
T = P~^, but also the coordinate representation by using 
the system energy E. Thus, the existence of first-order 
phase transition is not necessarily an intrinsic charac- 
teristic of a given system because of it depends on the 



external conditions which have been imposed. 

What happen when the ensemble inequivalence or any 
other anomaly is present at a given point p € 3 with- 
out matter the coordinate representation TZj e C used 
in the macroscopic description? Such anomalies are not 
representation-dependent, that is, they do not depend 
on the external conditions imposed to the system. These 
kind of anomalies are intrinsic characteristic of the sys- 
tem, and consequently, they have a topological character. 
A anomaly of this kind reflects the ocurrence of consid- 
erable changes in the microscopic picture of an isolate 
Hamiltonian system, i.e.: an ergodicity breaking J3Cll |. 

This idea brings us the recently proposed hypothesis 
about the topological origin of the phase-transitions [SJI ■ 
According to the ideas developed by these authors, the 
origin of the chaoticity of the microscopic dynamics and 
the origin of the phase transitions at the macroscopic 
level delay on the geometrical and topological structure 
of the configurational space. We think that such abrupt 
topological changes in the configurational space could 
be took place when they are associated with topologi- 
cal anomalies at the macrocopic level, that is, when such 
anomalies are not dependent of the coordinate represen- 
tation of the space 9. The non correspondence between 
top olog ical changes and phase transitions observed in 
ref. 36] might be explained due to the non topological 
character of such phase-transitions. The following theo- 
rem establishes the non existence of topological ensemble 
inequivalences in those systems controlling by only one 
scaling invariant parameter e, Ta (s) — e. 

Theorem 3 Let s (e) be the entropy per particle of a 
given system which is controlled by only one scaling in- 
variant variable e. Let us also suppose that the second 
derivative of s (e) exists, being s (e) a convex function 
in the interval Si < e < 62, and a concave function else- 
where. There is a reparametrization change ip : TZ^ — > TZip 
in which the entropy s becomes a concave function. 



function in the new coordinate representation TZ^p when- 
ever a (e) > 0. This demand is very easy to satisfy by 
considering: 

f \ _ ( d^s (e) /de^ wherever s (e) be convex 
^^' 1 —d'^s{e)/de^ wherever s (e) be concave ' 

(40) 
which directly leads to the following bijective map: 



de 




(41) 



where we have considered here the possibility of intro- 
ducing a different values for integration constant C in 
different regions. By choosing Ci = 1, C2 = l//3^ (ei) 
and C3 = /3^ (^2) //9^ (ei), we not only ensure the con- 
cavity of the entropy, but also the continuity of the first 
and second derivatives of the bijective map ^p{£). ■ 

The analysis about the existence of topological ensem- 
ble inequivalences for two or more controlling variables 
is an open problem. While topological ensemble inequiv- 
alence are not present in the unidimensional case, they 
are not the only topological anomalies presents in this 
context. The existence of a discontinuity in the first 
derivative of the entropy in the thermodynamic limit is 
a example of anomaly which could not be avoided by 
a reparametrization change Lp (e) with a continue first 
derivative dip (e) /de. Such anomalies have been recently 
reported in the astrophysical context, and they are usu- 
ally classified as microcanonical phase-transitions |37l| . 
Existence of singularities in the second derivatives of 
the entropy are also unavoidable by reparametrizations, 
which are related with the well-known second-order phase 
transitions appearing due to the existence of an sponta- 
neous symmetry breaking at the microscopic level 30j. 
Interestingly, the above topological anomalies appear due 
to an ergodicity breaking in underlying the microscopic 
dynamics [3(1 13^. 



Proof. The transformation rule of the second derivative 
of the entropy during the reparametrization change (p : 
TZ^ -^ TZ^p is given by: 



de \ d^s d^e ds 



dip^ \dip J de^ dip^ ds 
Considering the auxiliary function a (e): 

oip'^ \op J 



(37) 



(38) 



the first derivative of the bijective map ip (e) could be 
rewritten as follows: 



— Cp (ejexp ' ' 



de 



die) 



(39) 



where C is certain positive integration constant and 
P (e) = ds (e) /de. The entropy s will be a concave 



V. ASTROPHYSICAL SYSTEM 

Let us to illustrate the relevance of the scaling selfsim- 
ilarity and reparametrization invariance by considering 
the microcanonical description of a nontrivial nonexten- 
sive Hamiltonian system: the selfgravitating gas of iden- 
tical nonrelativistic point particles interacting through- 
out the Newtonian gravity: 



Hn =Kn + Vn 



-^ 2m^* 2 ^^ 



1 



('¥i> 



(42) 

The long-range singularity of the gravitational potential 
will be avoided by enclosing the system in a spherical 
rigid container of radio R. The short-range singularity 
will be regularized by considering a mean-field approx- 
imation. The above conditions leads to the well-known 
isothermal model of Antonov [Sa . 



After the integration of the hnear momentum variables 
p's, the accessibe phase-space volume can be expressed 
by: 



W ^Se- 



1 



2rm: 



Nl V h^ 



4jv 



j37V^ 



1 



r(iiv: 



{E-VnM) 



|Af-l 



where {r} = (ri, r2, . . . r„), (f^r = Hk'^^^k, and Se a 
very small energy constant. Introducing the character- 
istics units for the linear dimension Iq = R, total mass 
M = Nm, and total energy Eq = GM^ /R, and consid- 
ering a very large system size N in order to consider the 
Stirling approximation F (x) ~ (x/e) , the above expres- 
sion could be rewritten as follows: 



1 

iV 



■InM^: 



■ln/CoRN^\+lnuj{U), (43) 



where Co = ATrGm^/Sesh^ and uj (U) is given by: 

Inc. ([/) = Km 1 In ( fd'^^ {U - V^ U})^^ 

' (44) 

being U = E/Eq and (, — r/l^. The entropy per par- 
ticle (|43|) is scaling invariant under the following scaling 
selfsimilarity transformations: 



N{a)^aN,R{a)=a-^R^ E{a)=a^E, (45) 
result supporting the thermodynamic limit: 

N -^ oo, keeping fixed — - and RN^ (46) 

already obtained in refs.|28l l29l |. Now on so{U) = 
In a; (U) will be referred simply as the system entropy. 

The progresive calculation of (|44|l is carried out by con- 
sidering a mean-field approximation and them using the 
steepest descend method whose details will be omitted 
here in the sake of brevity. Such procedure leads to the 
following min-max problem: 

so (U) ~ minmax ( ^ In ([/ - 1/ [p]) + s [p] + fi {1 - N [p]) 

(47) 
where the functionals V [p] (potential energy) , s [p] (local 
entropy) and N [p] (particle number) are given by: 

v[p] = -lld^id^e^f^, N[p]^Jd^^p{o, 

^[P]-Id'^ {-piOlnpiO+pm- (48) 

The problem H47|l leads to the following Boltzmann 
distribution for the particles density p: 



p {^) = exp (-^ - P^ {^)) , 
where /? is the canonical parameter: 
dliiuj{U) 3 



P 



dU 



2{U~V)' 



(49) 



(50) 



/i, the chemical potential associated with the normaliza- 
tion condition: 



Jd^CpiO^h 



(51) 



and if (^) is the gravitational potential related with p 
throughout the Poisson problem: 



A^(0=47rp(0, ^(1) = -1, ^'(1) = 1. 



(52) 



The equation (|50|l can be rewritten in order to obtain the 
caloric curve U versus [3: 

U^^ + V,wheTeV^^Jd'^p{OfiO- (53) 

The solution with spherical symmetry of the above 
problem is numerically implemented by introducing the 
function $ (^) = —p + ln(3 — l3(p {£,), and solving the fol- 
lowing Poisson-Boltzmann problem: 






exp(<l>(0), (54) 



whose boundary conditions at the origin are: 



d 






(55) 



where V' — lii(/3pc), being p^ in the central density 
Pc = P (0). Thus, the functions p (^), Lp (^) as well as the 
macroscopic variables and thermodynamic potentials f7, 
/?, /i, So are obtained as functions of the parameter "0, 
i.e.: /3 and p are given by: 

$ (1; V) = -^ + ln/3 -f /3, ^<i> (1; V) = -p- (56) 

The FIG n shows the microcanonical thermostatistical 
description of the Antonov isothermal model. The caloric 
curve (panel a) and the central density Pc versus U de- 
pendence (panel b) show the existence of thermodynami- 
cal states with a negative heat capacity (branch A-B), as 
well as the ocurrence of a gravitational collapse pc —^ oo 
for low energies. The branch A-B-C corresponds to 
the equilibrium thermodynamic states, while the other 
branch are just unstable saddle points. According to the 
caloric curve, there are no equilibrium thermodynamical 
states when U <Ua^ -0.334 or /3 > /3b = 2.52. The 
ensemble inequivalence in the branch A-B is associated 
with the convexity of the entropy (panel c) in this region. 

As already discussed in our previous paper |l3, when 
this model system initially isolate with a total energy 
U £ {Ua,Ub) (branch A-B) is put in a thermal con- 
tact with a heat bath with /3 > /3b, it becomes unsta- 
ble developing an isothermal gravitational collapse [38| . 
so that, the Antonov isothermal model is very sensible 
the external infiuence of a heat bath in the branch A-B. 
Such instability can be associated with the existence of 
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FIG. 1: Themiostatistical description of the Antonov isother- 
mal model by using the representation TLij- 



FIG. 2: Thermostatistical description of the Antonov isother- 
mal model in the representation TZ^p. 



a first-order phase transition from a gaseous phase (with 
Pc finite) towards a collapsed one (pc -^ oo). 

The above thermodynamical description was per- 
formed in the TZu representation. According to the Theo- 
reni|31 the entropy convexity could be avoided by consid- 
ering an appropriate reparametrization change. In order 
to obtain such appropriate representation, let us to ob- 
serve the ^ versus U dependence also shown in panel d 
of the FIG^ It is easy to see that, dismissing the unsta- 
ble branch, there is a biunivocal correspondence between 
these macroscopic observables for all relevant equilibrium 
thermodynamical states. This fact supports that we can 
obtain an appropriate representation TZ^p by considering 
certain experimental arrangement which keeps fixed the 
observable x = exp (ip) — PPc- The map ip : TZjj -^ TZ^ 
can be established by taking into account that x is the 
canonical parameter in the new representation, which is 
related with (3 by the transformation rule (|28|l : 



X = exp [tp) = —P 



dip 



(3 [ip) exp (-'0) 



dUjtp) 

dip ■ 
(57) 



FIG 121 shows the thermostatistical description of the 
Antonov isothermal model in the representation TZi^,. The 
reparametrization change Lp : TZjj — > TZ^p obtained from 



the numerical integration of the equation (|57(l is illus- 
trated in the panel a (which is evidently a bijective map), 
while the "caloric curve" x versus Lp and the entropy 
So versus ip dependence are shown in the panel b. The 
reader can notice that no ensemble inequivalence takes 
place in the representation TZ^ , and consequently, the mi- 
crocanonical description of this model system in the ther- 
modynamic limit becomes equivalent to the one carried 
out by considering the following " canonical ensemble" : 



^B ix, N) 



Z{X,N) 



exp 



'XNp ( —Hn 



(58) 



associated with an experimental arrangement which 
keeps fixed x = f^Pc- Since no ensemble inequivalence 
is observed, no topological first-order phase transition is 
present in the Antonov isothermal model. 

It is very important to realize that the thermostatisti- 
cal description obtained from the numerical integration 
of the Poisson-Boltzmann problem (|54|l makes use (in an 
implicit way) of the above " canonical ensemble" since the 
spatial functions p and p, as well as the thermodynamic 
variables and potentials are obtained as functions of the 
parameter ip = lnx- 



VI. CONCLUSIONS 

We have shown in the present work that an appro- 
priate thermostatistical description of the nonextensive 
Haniiltonian systems could be performed by consider- 
ing their scaling selfsimilarity properties as well as the 
reparametrization invariance of the microcanonical de- 
scription. 

As already discussed in previous sections, the 
reparametrization invariance is just an internal symme- 
try of the microcanonical description appearing as a con- 
sequence of its dynamical origin. This symmetry allows 
us to satisfy the necessary conditions i) and ii) (see in 
section Ullll for the ensemble equivalence in the thermo- 
dynamic limit N —> oo, and the consequent well-defined 
character of the thermodynamic formalism based on the 
Legendre transformation H12I) among the thermodynamic 
potentials. The consideration of the reparametrization 
invariance makes possible the developing of a geometri- 
cal thermodynamic formalism in which the ordering in- 
formation based on the concavity of the Boltzmann en- 
tropy could be defined in a topological fashion, suggesting 



in this way a new vision of the phase transition concept 
which seems to be related with the called Topological Hy- 
pothesis about the topological origen of the phase tran- 
sitions |3J|. 

Our analysis of the Antonov isothermal model shows 
the relevance of the selfsimilarity and the reparametriza- 
tion invariance in the context of the astrophysical sys- 
tems. An interesting remark derived from this study is 
that the reparametrization invariance could be used as 
a powerful tool in order to extend a canonical-like de- 
scription (|58|l for those thermodynamical states with a 
negative heat capacity, fact which leads to a natural im- 
provement of the well-known montecarlo methods based 
on the consideration of the canonical weight exp {—f3E) 
|39| . A very important problem is still open: A method- 
ology to derive the relevant scaling selfsimilarity in a 
given application. Omitting some exceptional cases like 
the Antonov problem reconsidered in the present paper, 
derivation of the selfsimilarity scaling transformations is 
a very difficult task which demands an extension of the 
Renormalization Group methods [30j. 
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